Topological indices, defects and Majorana fermion in chiral superconductors 
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We study theoretically the role of topological invariants to protect the Majorana fermions in a 
model of two-dimensional chiral superconductors which belong to class D of topological periodic 
table. A rich phase diagram is revealed. Each phase is characterized by the topological invariants 
for 2d (Z) and Id (Z2), which lead to the Majorana fermion at the edge of the dislocation and the 
core of the vortex. Interference of the Majorana fermions originated from the different topological 
invariants is studied. The stability of the Majorana fermion with respect to the interlayer coupling, 
i.e., in 3d, is also examined. 
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Topological classification of the electronic states in 
solids has shed a new light on the band structure of 
solids and also the superconductivity [l[. Initiated by 
the proposal of Z2 topological invariant and the quantum 
spin Hall effect in two-dimensional spin-orbit coupled sys- 
tem [2] and its extension to three dimensions 3] , the more 
generic topological classification scheme, i.e., topological 
periodic table, is now available based on the three sym- 
metries, i.e., time-reversal (0), particle-hole (H), and chi- 
ral (II) symmetries [4j, [fj| . Here S symmetry is due to the 
superconductivity in the usual situation. There are 10 
classes, and the topological invariant is specified depend- 
ing on the dimensionality d of the system to character- 
ize the nontrivial topological states. Later, this "10-fold 
way" has been extended to the classification including 
the textures such as domain walls and dislocations [fj . In 
this case, the dimensionality D of the real space manifold 
surrounding the texture plays a key role, and 8 = d — D 
is replacing d for the topological classification [6] . In the 
topological periodic table, one can recognize the period- 
icity called Bott periodicity which relates the different 
classes in the "diagonal" direction. This periodicity can 
be understood by the continuous mapping of the Hamil- 
tonian connecting the different classes and different 8 by 
one. On the other hand, one can also consider the con- 
nection in the "horizontal" direction, i.e., dimensional 
reduction [7|. As an example, class All has been charac- 
terized by Z2 both for (5 = 3 and 8 = 2. The former cor- 
responds to the Z2 invariant corresponding to the strong 
topological insulator (TI) in 3d, while the latter to the 
quantum spin Hall system in 2d (D = 0) or the weak 
TI in 3d (D = 1). Namely, a 3d system is characterized 
by 4 Z2 topological invariants v§\v\v<iv<j n where vq = 1 
indicates the strong TI, while vq = with at least one 
nonzero ^1.2,3 means the weak TI. These ^1,2,3 are the 
topological invariants for (5 = 3 — 1 = 2, and guarantee 
the existence of gapless one-dimensional mode along the 
dislocation m. 



This topological periodic table provides the powerful 
guiding principle also for the topological superconductors 
(TSs) . Especially, the Majorana fermions expected to ap- 
pear at the edge or the core of the vortex in TSs attract 
intensive interests from the viewpoint of the quantum 
information technology |9Ml4 1 . Therefore, it is an impor- 
tant theoretical issue to design the Majorana fermions 
in realistic systems. Proximity-induced superconductiv- 
ity in 3d T I llal. the superconductivity in a doped TI 
Cua;Bi2Se3 [16lll7|. the possible TS in noncentrosymmet- 
ric systems with the Rashba spin splitting [18l423j , and p- 
wave superconductivity in Sr2Ru04 24j are the promis- 
ing candidates as the host of the Majorana fermions. As 
pointed out in ref. [25j , most of the theoretical proposals 
for the Majorana fermions are based on the two mod- 
els, i.e., p-wave pairing in the one-dimensional spinless 
fermions (Kitaev model [9() and the p + ip pairing super- 
conductor. 

In this paper, we study theoretically the topological 
invariants and their relation to the protected Majorana 
bound states in a model of class D chiral superconductors 
containing both the Kitaev model and p + ip supercon- 
ductor in the limiting cases. The topological invariant 
for class D is for (5 = 3, Z for <5 = 2, and Z2 for 
(5 = 1. Therefore, there is no "strong TS" in 3d, while 
the 2d system is characterized by Z topological invari- 
ant and the Id system by Z2 topological invariant. The 
purpose of the present Letter is to reveal the topologi- 
cal phase diagram characterized by these invariants, and 
the associated Majorana fermions at the textures such as 
dislocations and vortices. 

We consider a generalized model of p + ip wave super- 
conductor on a square lattice in 2d. The Hamiltonian 
can be written as H = J2k^k^(^)^ k with 



H{k) 



2t x cos k x -\-2t y cos k y —fi d x sin k x —id y sin k y 
d x sin k x -\-id y sin k y fi—2t x cos k x —2t y cos k y 



(i) 

and C£ = (cj^c-k)- This 2x2 Hamiltonian matrix can 



be expressed as H(k) = H(k x ,k y ) — h(k) ■ cr where 
<j = (a x ,(j y , <r z ) is the vector of Pauli matrices. Since 
(C*l k ) T = a x C k , H(k) should satisfy 



T„x 



H(k) = -a x H{-kYa 



(2) 



where T means the transpose. This condition leads to 
the relation |25| 

h x ,y{k) = -h XtV (-k), h z (k) = h z (-k). (3) 

Therefore, for the time-reversal invariant momenta 
(TRIM), which satisfy k = — k, only h z (k) cannot be 
nonzero, i.e., h{k) points either in +z or — z directions 
as long as the gap opens, i.e., \h(k)\ > 0. There are 4 
TRIMs in this 2d model, i.e., k a = (0,0), (tt,0), (tt,tt) 
and (7r, 7r), and the sign s a = ±1 of the corresponding 
h z . As will be discussed, s a determines the Z 2 topologi- 
cal invariants and the parity of the Z invariant. 

Now let us start with the Z 2 invariant. For this 
purpose, let us consider the Id Hamiltonian with fixed 
k x = ir in cq.([T]), i.e., 



H(k ns =TT,k y ) — 



-2t^ J r2t y cos(k y )— fi — idy sin(fcy) 

id y sin(ky) fi-\-2t x — 1t y cos(k y ) 



(4) 

which is nothing but the Kitaev model for a one dimen- 
sional topological superconductor [9j . The ^^topological 
invariant v x is related to the "polarization" [7(. 



-£=P(k x ) = / -—^-a y (k x ,ky) modi 

1 J —K ^7T 



(5) 



given by the Berry phase vector potential aj(k x ,k y ) = 

H) T, n -.occu^diif\ d / dk i\ nk )^ and is s iven b y (- 1 )"* = 

S( 7 r i o)S(7r .-n) [25l |26|. Therefore, we can easily obtain the 
Z 2 topological invariant v x as 



1 for % + %\ < 1^1 
for \t x + $\ > \t y \ 



(6) 



The topological invariant v' x for k x = can be also cal- 
culated in a similar way. From these equations, it is clear 
that strengths of d x and d y are not related to the topolog- 
ical numbers if they have finite values. The Z 2 invariants 
v y , v' are obtained in the similar way. 

On the other hand, the Z topological invariant v is 
nothing but the Chern number, i.e., the wrapping num- 
ber of the mapping from the 1st Brillouine zone of k to 
the unit sphere h(k)/\h(k)\ and is given by 

u= — 7, — -{9k x a v (k x ,ky) -dk y a x (k x ,k y )} . (7) 

JJbz ^ 

Equations © and © lead to the relation @, E0, EH, El] 



t x + L = - 




*x ~r t. 



FIG. 1: The topological phase diagram of a model in eq. |T]) 
for chiral superconductors in 2d characterized by v : u x u y . 
The lines where the energy gap closes divide the (i^i^-plane 
into nine domains. The system is the strong topological su- 
perconductors (y — 1) in the domains II, III, VII and VIII, 
while it is weak topological superconductor in the domains I, 
IV, VI and IX with v — but some of Z2 invariants being 
nonzero. In the domain V, the system is the trivial strong 
coupling superconductor. 



In summary, our model is characterized by Z topolog- 
ical invariant v and two Z 2 topological invariants v x and 
v y . This is the general result, and superconductors in 
class D in 2d are characterized by v : v x v y . From the 
above consideration, these topological invariants depend 
on the hopping integrals t x , t y and the chemical potential 
fj,, while do not depend on the pairing amplitudes d x , d y 
as long as they are finite. Therefore, we show in Fig. [1] 
the phase diagram of the present model in the plane of 
(t x , t y ) for fixed [i. The lines where the energy gap closes 
divide the (t x , t y ) -plane into nine domains. Electronic 
states are characterized by topological invariants in each 



v mod 2. 



(8) 



domain. When ^ is larger than \t x \ + \t y \, i.e., domain 
V, the pairing state is topologically trivial since it cor- 
responds to the strong coupling limit. In the domain 
II, IE, VH and VHI, electronic states have both the Z 
and Z 2 topological invariants. Note that the sign of the 
chemical potential [i matters for the Z 2 invariants. In 
the domain I, IV , VI and DC, the anisotropy between t x 
and t y is large and hence the system behaves basically 
as the weakly coupled chains of Id Kitaev models, and 
the system has only Id Z 2 topological invariants but Z 
topological invariant v = 0, so they are weak topological 
states. 

Now we turn to the consequences of the topological in- 
variants. The Z 2 invariant v x (u y ) ensures that the prop- 
agating Majorana fermion channels appear at the edge 
along x direction (y direction). They have zero-energy 
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FIG. 2: The edge dislocations (indicated by yellow points) 
in two dimensional system and three dimensional system. 3d 
system is made by stacking 2d systems and adding a hop- 
ping integral along the z-direction with the edge dislocation 
isolated on a layer. 



states at k x = it (k y — tt). Also the edge dislocation of- 
fers the stage for the zero energy Majorana bound state 
when the following equation is satisfied (see Fig. [5J. 



BG = 1 



mod 2. 



(9) 



In this equation, we define G = ■^-{v x h x J \-v y h y ). h x and 
b y are reciprocal lattice vectors of x- and y-direction, re- 
spectively, and B is the Burgers vector characterizing the 
dislocation. We numerically calculate these zero-energy 
states as follows. To introduce a periodic boundary con- 
dition, we introduce two edge dislocations with the Burg- 
ers vector B = ±e x . We represent two edge dislocations 
by adding lattice sites between them. Edge dislocations 
are separated by a half system size. Calculations were 
done on a 40 x 40 unit cell system with a periodic bound- 
ary condition along x- and y-directions. The parameters 
are t x = 0.5, t y — 0.5, d x — 0.6, d y = 0.6, fi — —0.2. In 
these parameters, the topological invariants are 1 : 11. 
In this case: G = ^(t>x + bj,) and B = ±e x , B ■ G = 1 
is satisfied, so zero-energy states appear at edge disloca- 
tions. 

The results of our numerical calculations are shown 
in Fig. [3J The left figure of Fig. [3] indicates the energy 
levels of our model in the presence of edge dislocations. 
It is clear that zero-energy states exist, which are 2-fold 
degenerate because the 2 edge dislocations are present 
in the system. The right figure of Fig. [3] indicates the 
probability distribution of the zero-energy states. Zero- 
energy states are localized at each edge dislocation. 

In the presence of edge dislocations with the Burgers 
vector B = ±e x , zero-energy states appear in the do- 
main I,n,VH and IX in Fig. Q] in the case of /x > and 
zero-energy states appear in the domain I, m, VIII and 
IX in Fig. [1] in the case of /J, < 0. In particular, zero- 
energy states in the domain / and IX can be intuitively 
interpreted as follows. The weak topological supercon- 
ductors in these domains are adiabatically connected to 
a stack of the Kitaev models [9( for a Id topological su- 
perconductor along y-direction. In the presence of edge 



FIG. 3: The energy levels and the probability distribution of 
zero-energy states at each site in the presence of edge dis- 
locations: The left figure indicates the energy levels. The 
right figure indicates the probability distribution for the zero- 
energy states. The zero energy Majorana bound state appear 
at each of the edge dislocations. 



dislocations, the edges of Id topological superconductor 
appear at edge dislocations as shown in Fig. [2] so zero- 
energy states appear there. In general, the existence of 
zero-energy states is proved by the same method as in 
ref. [1. 

Next we consider the interference of the Z and Z 2 topo- 
logical invariants, which is realized in the present model. 
When Z is nonzero, the zero energy Ma jora na bound 
state is realized at the core of the vortex [12|, ll_3( . It is 
expected that, if the dislocations are in the crystal, they 
act as the pinning centers of the vortex, and hence there 
are two reasons for the existence of the Majorana bound 
states when Z 2 invariant is 1. This situation occurs in 
the domains II, HI, VII, and VIE in Fig. 1, and the in- 
terference of these two mechanisms is an issue. Figure 0] 
summarizes the calculated results, in which there are dis- 
locations and vortex cores at same positions. The proba- 
bility distributions are plotted for the zero energy states if 
any. To introduce a periodic boundary condition, we in- 
troduce two vortices with the winding number 1 and two 
vortices with the winding number -1. We consider the 
two cases of topological invariants 1 : 00 (upper panels) 
and 1 : 11 (lower panels). The parameters are t x = 0.5, 
t y = 0.5, d x — 0.6, d y — 0.6, /1 = 0.2 for the former, 
while they are t x — 0.5, t y — 0.5, d x = 0.6, d y = 0.6, 
/j = —0.2 for the latter. Calculations were done on a 
40x40 unit cell system with a periodic boundary con- 
dition along x and y-directions. There are dislocations 
and vortex cores at same positions in the right panels 
while only dislocations are there in the left panels. They 
are separated from each other by a half system size. In 
the case of the topological invariants 1 : 00, zero-energy 
states do not appear at edge dislocations. Zero-energy 
states appear when the dislocations and vortices exist at 
the same time because of the Z invariant v — 1 and vor- 
tices. We have also confirmed that the zero-energy states 
appear with only the vortices. In the case of the topo- 
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FIG. 4: The interaction between zero-energy states localized 
at edge dislocations and zero-energy states localized at vor- 
tex cores: This figure shows the probability distribution of 
zero-energy states at each lattice sites in each case. The 
interaction between zero-energy states at edge dislocations 
and zero-energy states at vortex cores eliminates zero energy 
states when dislocations and vortex cores exist at the same 
positions. 



logical invariants 1:11, zero-energy states appear when 
dislocations exist but vortices do not exist. The interac- 
tion between zero-energy states at edge dislocations and 
at vortex cores eliminates zero energy states when they 
coexist at the same position. This can be naturally un- 
derstood that the two Majorana bound states due to Z 
and Z2 invariants interact with each other as approach- 
ing to each other, and lift the degeneracy to have finite 
energies. 

We generalize this model into the model in 3d. The 
3d model is constructed as a stack of 2d models with 
a hopping integral t z along the z-direction as shown in 
Fig. 2. t z causes the finite region where the energy gap 
closes in the (t x ,t y ) -plane. The energy gap closes in the 
domain which satisfies \t x +t y — ^\ < \t z \, \— t x +t y — f | < 

The 



II < 
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-t. 



' u 



£ < t z 



l^zl) {"x ty 2 I ^ I L z I ^ L L x L y ■> 

domain where the energy gap opens is connected to the 
domain of t z = 0, so electronic states in the domain have 
the same topological invariants as the 2d model, while 
the topology is trivial for class D in 3d. We introduce 
the similar defect as the 2d system as indicated in Fig. [3] 
These defects are not edge dislocations. One layer has 
the same defect as the 2d system and other layers do 
not have it. We ensure that zero-energy states appear at 
those defects in the same manner as the 2d system by 
the numerical calculations. The parameters are t x = 0.5, 
t y = 0.5, t z = 0.1, d x = 0.6, d y = 0.6, ju = -0.5. In 
this case, the electronic state is adiabatically connected 
to the state in the domain 77/ in the case of fx < 
in Fig. [TJ If t z = 0, zero-energy states appear in these 
parameters because the state in the domain 77/ in the 
case of /i < in Fig. [1] has topological invariants 1 : 



11 and one layer has edge dislocations with the Burgers 
vector B = ±e x . Calculations were done on a 10x20x10 
unit cell system with a periodic boundary condition along 
x, y and z-directions. We found that these zero-energy 
states survive a finite- value z direction hopping integral t z 
unless the energy-gap closes. This edge dislocation is Od 
object, and we need D = 2 sphere to enclose this defect in 
3d, and hence (5 = 3 — 2 = 1. Therefore, the Z2 invariant 
protects the existence of the Majorana zero energy bound 
state. This is an example of the "weak- weak" topological 
superconductor, where the topological invariant reduced 
by 2 dimensions is relevant. This indicates that Majorana 
fermions can exist even in a 3d system. 

Now we discuss the relevance of the present results to 
the real materials. There are several candidates for the 
chiral superconductors. The 2d Rashba system in a semi- 
conductor with the proximity to the s-wave superconduc- 
tor and the ferromagnet is a promi sing can didate, which 



imising 

MM 



shows the spinless p+ip pairing 19, 20, 22, 23]. The elec- 
tron density is considered to be small and concentrated 
near the r-point of the 1st Brillouine zone. Therefore, 
usually the continuum approximation is used to describe 
this system. In our phase diagram (Fig. [T]) , this situa- 
tion corresponds to the domain VII with t x = t y < 
with negative chemical potential, which is characterized 
by v — 1 and v x = v y = 0. Therefore, we do expect the 
Majorana bound state at the core of the vortex, while 
there is no Majorana bound state at the edge disloca- 
tion. To realize more interesting situation where both 



and 



are nonzero, it is required to have the Fermi 



pocket near the zone boundary or zone corner. This is, in 
principle, possible once the sign of the hopping integral is 
reversed. Another interesting candidate is Sr 2 Ru04 [2J], 
which is believed to be a quasi-2d p + ip superconduc- 
tor with the non-zero Z invariant. In this system, there 
remains the spin degeneracy of the Fermi surface, and 
hence the Hamiltonian matrix is at least 4x4 instead of 
2x2 discussed in this paper. This means that the Zi 
invariants are zero due to this degeneracy, although Z 
invariant can be nonzero. 

To summarize, we have studied a model of chiral su- 
perconductors including both the Kitaev model in Id and 
p + ip superconductor in 2d as limiting cases. This model 
shows a rich phase diagram (Fig. 1) characterized by the 
Z and Z 2 topological invariants, which control the ap- 
pearance of Majorana bound states at the edge disloca- 
tions and vortex cores. This offers an explicit case where 
the topological periodic table can be successfully applied 
including the topological textures, and the presence of 
the Majorana bound states are shown also numerically. 
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